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Abstract. 1/f noise, the major source of dephasing in Josephson qubits, may be 
produced by an ensemble of two-level systems. Depending on the statistical properties 
of their distribution, the noise distribution can be Gaussian or non-Gaussian. The 
latter situation is realized, for instance, when the distribution of coupling strengths 
has a slowly decaying power-law tail. In this regime questions of self-averaging and 
sample-to-sample fluctuations become crucial. We study the dephasing process for a 
class of distribution functions and analyze the self-averaging properties of the results. 
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1. Introduction 

In Josephson qubits dephasing is dominated by low- frequency noise, often with a 1/f 
power spectrum, due to fluctuations of background charges, magnetic fluxes, or critical 
currents pfl El Oj- While irrelevant for the relaxation process with time scale 7\, low- 
frequency noise dominates the dephasing time T 2 *. Standard NMR echo techniques 
allow one to reduce dephasing by rendering the low frequency spectrum ineffective [JJ. 
Operation at optimal bias points, chosen such that the linear longitudinal coupling of 
the qubit to the 1/f noise source vanishes, proved to be very successful in increasing 
the dephasing time [2]. Further progress in this direction may require an improved 
understanding of the mechanisms causing 1/f noise and of its statistical properties. It 
was realized recently that qubits themselves can be used to study the noise properties 
of their environment jH E] , and an interesting relation between the low-frequency 1/f 
and the high-frequency charge noise was observed 6J. An extensive study of dephasing 
due to both charge and flux noise was undertaken in Ref. [Zj. 

Still, many questions remain open. If the number of fluctuators contributing to 
the 1/f noise is large, one could expect Gaussian statistics [JJ IE]- In Ref. jHj and 
following work the role of individual, strongly coupled fluctuators was emphasized, and 
it was suggested [TU] that even ensembles of many fluctuators may produce strong 
non-Gaussian effects, emerging as a result of rare configurations in which dephasing 
is dominated by a small number of very strongly coupled fluctuators. As far as we can 
judge, the decay laws observed in Ref. [Zj cannot be fully explained by either of these 
theories. 

As the experiments are performed on individual systems with a particular 
configuration of the fluctuators, it is important to understand whether the predicted 
decay laws are self-averaging or have strong sample-to-sample fluctuations. Here we will 
analyze a class of distribution functions for the coupling strengths of the fluctuators. 
We determine the ensemble-averaged decay laws (extending the results of Ref. [TU] ) 
and analyze which of them are self-averaging. We study both dephasing due to 
linear longitudinal coupling and dephasing at the optimal point where the coupling 
is quadratic. 

2. 1/f noise from two-level fluctuators (TLF) 

1/f noise is often attributed to a collection of bistable systems, switching randomly 
between two states [TTJ. On one hand, such a model provides a natural explanation of 
1/f noise. On the other hand, in many samples distinct two-level fluctuators (TLF's) 
were detected. In metals this switching causes conductance fluctuations O H3j and, 
consequently, 1/f noise of the transport current. In Josephson junctions it causes the 
critical current to fluctuate f3J E3- More generally, spin bath environments were 
analyzed in Ref. JUJ. In charge qubits the TLF's contribute to the fluctuations of 
the gate charge controlling the qubit. The TLF's are characterized by their coupling 



3 



SnH = ^2 • ( 2 ) 



strengths to the qubit, v n , which may vary depending on the location of the respective 
TLF. The fluctuating quantity that couples to the qubit, X(t), contains contributions 
from all TLF's: 

X(t) = ^V n <T rhz (t) (1) 
n 

Each fluctuator switches randomly between two positions, denoted by a n ^ z — ±1, with 
rate 7„ (for simplicity, we assume equal rates in both directions for the relevant TLF's) 
and thus contributes to the noise power Sx{oj) = Yin Sn(u) with 

A set of TLF's produces 1/f noise when the switching rate 7 depends exponentially 
on a physical quantity, I, with a smooth distribution. For instance, 7 oc e~ l ^ l °, 
with I distributed uniformly over a range much wider than l Q , translates in a log- 
uniform distribution of the switching rates, with probability density Piff) oc I/7 in 
the corresponding exponentially wide range 7 min C 7 <C 7max- In this range the total 
noise power thus scales as 

Sx{uj) oc / — — - oc — . 3 

J 7 Y + u 2 \uj\ 

An example is a particle trapped in a double-well potential, whose tunneling rate 
through the potential barrier depends exponentially on both the height and the width 
of the barrier, leading to 1/f noise. Another example is thermally activated tunneling 
with rate ^§e~ E l k ^ J , where E denotes an activation energy. In this way the 1/f power 
spectrum observed in metals can be attributed to a broad (much wider than k^T) 
distribution of activation energies [T7j . 



3. Distribution of coupling strengths, self-averaging 

The analysis of decoherence in the presence of many fluctuators requires the study of 
probability distributions of coupling strengths and switching rates. In each particular 
sample one deals with specific fluctuators, i.e., with a realization of the set of parameters 
v and 7, drawn from this distribution. One should distinguish between quantities 
averaged over a statistical ensemble of samples and the results for a specific sample. 
This difference is essential, if the quantity under consideration is not self-averaging, i.e., 
if it has considerable sample-to-sample fluctuations. Such a situation arises if a quantity 
is dominated by contributions from a small number of TLF's. 

In Ref . [TU] , a continuous distribution of the parameters v n and 7 n was considered, 
with a long tail of the distribution of coupling strengths v n , such that rare configurations 
with very large v n dominate certain ensemble properties. It arises, e.g., from a 
uniform spatial distribution of fluctuators on a d- dimensional surface and a power- 
law TLF-qubit coupling [TU], v(r) oc l/r b . This results in a distribution of coupling 
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strengths P(v) oc l/v l+d l b . The joint distribution P(v,j), defined in the domain 
[v min , oo] x [y m in, Imax] and normalized to describe N fluctuators is thus 

P(v,l) = -%- (4) 

Here /i = d/b > 0, c = 1/ ln(7 max /7 min ), and r\ = VmmN 1 '^. One can also allow for 
fluctuations of A, but this does not change the results significantly. 

We consider a d-dimensional volume of typical size r max around the qubit containing 
a uniform distribution of TLF's. The typical distance between the strongest (closest) 
fluctuator and the qubit thus scales clS ^min r ^ J 

{V/Nf' d ~ rw/iV 1 ^. On the other hand, 
since the coupling strength was assumed to decay as v (r) oc l/r b the relation between 
the strongest and weakest coupling strength is given by v maiX /v min = (r max /r min ) 6 . 
Combining both results we find that the typical maximal coupling strength scales as 
■^max ~ ^min^ 1 ^- This does not exclude the existence of fluctuators with v 3> in 
certain realizations, as the long tail of the distribution function suggests. 

As examples of averaging over the distribution of coupling strengths and switching 
rates we calculate the noise produced by the ensemble of fluctuators, 

S x (u)= / dt,d 7 P(t, l7 )_-TL (5) 
J Y + u z 

distinguishing two cases: In one case, for fi < 2 the integral over v diverges at the upper 

limit. Hence the noise is dominated by the strongest fluctuator (s). Thus the result is 

sensitive to the properties of one or a few fluctuators and is therefore not self-averaging. 

Estimates below are based on cutting the integral at v = 77 = but one has to 

remember that for a comparison with experiments averaging (including the averaging 

over 7) makes little sense, since only a few TLF's contribute. 

In contrast, for fi > 2 the integral is dominated by fluctuators with v < rj. The 
weak fluctuators are most important, and due to their large number the noise is given 
by a sum of many comparable independent contributions. Consequently, the result is 
self-averaging, i.e., in different samples or runs of the experiment with \i > 2 one should 
observe the same noise amplitude. 

We now summarize the typical/average results for the noise, retaining only the 
leading contributions 

2-kA f c 77 2 : n < 2 (typically) 

S x (uj) = with A = { M (6) 

M cN(v 2 ) : fJl >2 



V )v 
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For ji > 2 we defined the average coupling strength of the TLFs, (v 2 ) — jr J dvP( 
Note that © is only valid for frequencies 7 min <C |u;| <ti 7m ax- At lower frequencies, 
M < 7min, Sx(w) tends to a constant, whereas at higher frequencies, \u\ > 7 max , Sx(w) 
crosses over to a faster power law decay oc 1/u 2 . 
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4. Longitudinal and transverse noise coupling 

We consider a qubit controlled (for simplicity) by a single parameter A and Hamiltonian 
H qb = - l -H,{X)a. (7) 

After an initial preparation in a coherent superposition of the qubit 's eigenstates, 
the effective spin precesses under the influence of the static field Hq, set by the 
control parameter Ao- Coupling to the environment disturbs this evolution, leading 
to decoherence. In many cases the effect of the environment can be modeled by classical 
and quantum fluctuations of A(i) = Xo + X(t), where X(t) fluctuates. For instance, in a 
charge qubit electromagnetic fluctuations of the control circuit as well as the background 
charge noise influence the gate voltage which controls the qubit. 

To proceed we expand the Hamiltonian to second order in the perturbation X, 



H *=~2 



H(\) + dH °X+ d2H ° X \ 
^o(A ) + ^ r X + ^ — + 



a. (8) 



Introducing the notations D\ = (1/h) dH /dX and D\ 2 = (l/h) d 2 H /dX 2 , we find in 
the eigenbasis of H (X ) a: 

H qb = --h (u 01 a z + 5u z a z + 5uj ± a ± ) (9) 

where fuj 01 = \H (\ )\, 8uj z = D Xz X + D\ 2)Z X 2 /2 + and 5u± = D X< ±X + .... Here a± 
denotes the transverse spin components (i.e., a x or a y ). The coefficients D are related 
to the derivatives of u; i(A): 

d\ 

and 



Dx, z , (10) 



Thus, in general, the coupling of noise to the qubit contains both transverse (Su z ) 
and longitudinal (5uj±) parts, and both may have linear as well as higher order (e.g., 
quadratic) contributions. 

5. Bloch-Redfield theory 

For weak, short-correlated noise the dynamics of the two-level systems (spins, qubits) 
can be summarized by Bloch's equations fHl IS] i n terms of two rates: the longitudinal 
relaxation (depolarization) rate r x = T-j" , and the transverse relaxation (dephasing) 
rate T 2 = T^ 1 . Evaluated perturbatively, using the golden rule, the rates are given by 

Ti = 2 S &»±( W = w oi) = ^D 2 x ±S x (uj = lu q1 ) (12) 

and 

r 2 = ^r 1 + r v , (13) 
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where 

T v = i£*Ju; = 0) = \d\ z S x (uj = 0) . (14) 

The dephasing process (|13|) is a combination of depolarization effects (I\) and of the 
so called 'pure' dephasing, characterized by the rate T v = T 2 * _1 . The pure dephasing is 
usually associated with the inhomogeneous level broadening in ensembles of spins, but 
occurs also for a single spin due to the longitudinal low-frequency noise. 



6. Pure dephasing for Gaussian noise, [i > 2 

If there are sufficiently many fluctuators in the environment, the central limit theorem 
(CLT) applies, and the noise is Gaussian. More specifically, since the central limit 
theorem applies to a large collection of equally distributed random quantities, one needs 
to have a large number of TLF's of each (relevant) "kind" (i.e., for each pair v, 7). 
This implies a regular distribution of coupling strengths, so that the relevant physical 
quantities are not dominated by a few TLF's at a boundary of the distribution. In 
particular, the distribution (J3J) gives rise to Gaussian noise if \i > 2. We will discuss now 
the pure dephasing derived from such Gaussian noise. The random phase accumulated 
at time t, 



t 

A<f> = D x>z J dt'X{t') 





is then also Gaussian-distributed. Hence the decay law, due to longitudinal noise 
(coupling to a z ) in a free induction decay (Ramsey signal) is given by fnit) = 
(exp(iA0)) = exp(— (l/2)(A</> 2 )). Averaging here is over the different trajectories of 
X{t) in repeated runs of the dephasing experiment. We obtain 

t 2 „o f + °° duo „ , , „ ut 



f R (t) = exp 



~2 D W 2^ Msinc2 T 



(15) 



where sine x = sin x/x. If most of the noise power is concentrated at frequencies u 1/t 
(static noise), then one can approximate sine |«1 and obtain 



/**(*) = exp 



t 2 

n 2 ^ 2 



(16) 



where a\ = f_°°(duj/27r)S\(uj) is the dispersion of X. In general, for static noise with 
(not necessarily Gaussian) distribution function P(X) the Ramsey decay is given by 



f^\t) = j d(X)P(X)e tU ^ xt , (17) 

i.e., by the Fourier transform of P(X). Static noise corresponds to a situation when X 
is constant during each run of the experiment but fluctuates between different runs. 
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In an echo experiment, the phase acquired is the difference between the two free 
evolution periods: 

t/2 t 



+ A</> 2 = -D XiZ J dt'X(t') + D XtZ J dt'X(t') 



t/2 



which after averaging over the trajectories of X(t) gives 



r r * 2 n 2 f +C ° duJ a ( \ ■ 2 . 2 ut 

f E (t)=exp --D Xz —S x (u)sm -sine — 



(19) 



1/f spectrum: Here and below we assume that the 1/f law extends over a wide 
range of frequencies, limited by infrared and ultraviolet cut-offs, 

In A A 

S\(uj) = - — - = — , for uj h <C |co>| <C uj c . (20) 
M rl 

The infra-red cutoff uji r is usually determined by the measurement protocol, as discussed 
further below. The decay rates typically depend only logarithmically on uj- IT) and details 
of the noise power below uj ir are irrelevant to logarithmic accuracy. For most of our 
analysis, the same applies to the ultra-violet cut-off uj c . However, for some specific 
questions considered below, frequency integrals may be dominated by uj ~ uj c , and thus 
the detailed behavior near and above u c (i.e. the "shape" of the cut-off) is relevant. We 
will refer to an abrupt suppression above u> c (S(w) oc Q{uj c — |u;|)) as a 'sharp cut-off', 
and to a crossover at u ~ uo c to a faster decay 1/ui — > 1/u 2 (motivated by modeling of 
the noise via a set of bistable fluctuators, see below), as a 'soft cut-off'. 

For 1/f noise, at times t l/uj ir , the free induction (Ramsey) decay is dominated 
by the frequencies u < 1/t, i.e., by the quasistatic contribution [20], and (fTHjl reduces 
to 



f R (t) = exp 



-f 2 Dl z A(ln— + 0(r 

UJj r Z 



(21) 



Here the logarithmically large part of the exponent originates from a static contribution 
of frequencies to < 1/t. Indeed it can be obtained from Eq. (fTB|) with o\ = 
2 ^^(du /2-n)Sx(yj) = A ln(l/c<jj r t). This contribution dominates the decay of fait). 
For the echo decay we obtain 

f E {t) = exp [-t 2 D\ z A -In 2] . (22) 

The echo method thus increases the decay time only by a logarithmic factor. This low 
efficiency of the echo has its origin in the high-frequency tail of the 1/f noise, which, as 
we note, influences the results strongly. For 1/f noise with a low cut-off u c the integral 
in Eq. ()19j) over the interval uj < uo c is dominated by the upper limit. For instance, in 
the case of a sharp cutoff, i.e., S = (A/\u\)9(u c — uj), we obtain 

f E (t)=exp(-±Dl z Auj 2 t^ . (23) 

On the other hand, for a soft cut-off, which we expect when the noise is produced by 
a collection of bistable fluctuators with Lorentzian spectrum, the integral in Eq. (|19|) is 
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dominated by frequencies u c < uj < 1/t, and we find \ia f E (t) oc D\ z Au c t i . In either 
case, one finds that the decay is slower by a factor ~ {u c t) 2 or u c t, respectively, than 
for 1/ / noise with a high cutoff, tu c > D\ tZ A 1 / 2 . 



7. Individual fluctuators 



We consider a single fluctuator coupled longitudinally to the qubit, whose contribution 
to the level splitting, v n (t) = v n a ntZ (t), switches between ±v n . For this case the free 
induction (Ramsey) and echo decays have been evaluated in Refs. [OJEI]. In the limit of 
high effective temperature, i.e., when the transition rates in both directions are equal, 
the decay functions, obtained by averaging over the switching history of <7„ j2 (t), are 
given by 



f R ,n{ t ) = ^ nt (cos fJ. n t+— sin fi n t] , (24) 



and 



f E ,n(t) = e " 7 "* ( 1 + — sin /^ + nr(l - cos M ) , (25) 



where fi n = a/ (Dv n ) 2 — 7^ and D = D\ z . In order to derive these expressions we 

t 

introduce the averaged phase factors x±if) — ( ex P(*/ dt' D v n (t'))) , averaged over the 
switching histories ending at v n (t) = +v n or —v n , respectively. Their dynamics is 
governed by the rate equations 

X+ = iDv n x+ - lnX+ + InX- , 

X- = ~iDv n x- - InX- + lnX+ ■ (26) 

The solution for fn >n (t) = X+if) + X-{t) is obtained by solving the coupled equations 
for the initial conditions x± = 1/2, which gives Eq. ()24j) . Similarly, for more general 
protocols, we have to analyze phase factors (exp(i J* dt'D g{t')v n {t'))) with appropriate 
time dependence of g(t). In this case the first terms on the right hand side of Eqs. (|2l)|) 
are modified accordingly. For the echo experiment we obtain in this way Eq. (J2HJ). 

The decay produced by a number of fluctuators is the product of the individual 
contributions, i.e., fnit) = Il n /i?, n (t) and /s(^) = n n fE,n(t). If the noise is dominated 
by a few fluctuators (this includes the case of many fluctuators in total, but a few of 
them with similar rates 7) the fluctuations of X(t) may be strongly non-Gaussian. 



8. Non-Gaussian effects, < 2 

Since we consider uncorrelated TLF's the total decay of coherence is the product of all 
single-TLF contributions, f(t) = Y\ n f n if)) where f n {t) is given by (jUj) and (J23) for 
the free induction and echo experiment, respectively. In Ref. ^Hj an ensemble-averaged 
value of ln/(t), denoted as (ln/(t)) F , was calculated for /z = 1. Here (■■ -)f denotes the 
average over the distribution of coupling strengths and switching rates Both free 
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induction decay and the "phase memory decay" (a protocol similar but not equivalent 
to the spin echo decay) were analyzed in the regimes t < 7"^ and t > 7"^. Below we 
will generalize these results to the range < /i < 2. 

As discussed above, the quantity (In f(t))p is relevant for experiments with specific 
samples only if the sample-to-sample fluctuations of In f(t) are weak, i.e., if In f(t) 
is self-averaging. Then, the experimentally observable decay law f(t) would be well 
approximated by exp((ln f(t))p). In Ref. ^U] the self-averaging was numerically 
confirmed for the phase memory decay in the regime t < 7mL- Here we analyze the 
self- averaging in four regimes: for the free induction and the echo cases, both in the 
limits t < 7~aa, and t > %^ ax . Specifically, we evaluate the ensemble-average (In f(t))p, 
given by an integral over the (v , 7)-space. In some cases this integral is dominated by 
a range in the 'bulk' of the distribution, which contains many fluctuators on average; 
this indicates that sample-to-sample fluctuations are weak. In other cases, the integral 
is dominated by the boundary of the distribution, indicating that the studied quantity 
is not self-averaging. Our analysis confirms the conclusion of Ref. [10] . obtained in one 
regime: for the echo decay at short times t < 7"^. We show further that in all other 
three regimes investigated the dephasing law is not self-averaging. In the calculations 
we assume that v m i n , 7 min are very low frequency scales, and I ft always exceeds them. 

8.1. Free induction decay 

For short times, t < 7~„x! we are effectively in the static regime, and the ensemble- 
averaged free induction decay is described by 

(ln|/ fl (t)|) F oc-(D A ^^ . (27) 

This result is dominated by the fluctuators with strength of order v ~ v^ x ~ r\ and thus 
is not self- averaging. For an experiment with a specific sample the results should be fitted 
by a contribution of one ([2*1]) or a few fluctuators, rather than by the ensemble-averaged 
behavior ([27]) . We can also estimate the typical decay law for short times t < rf~ x . In 
every realization there will be a few strongest fluctuators, typically with strength f max - 
For t <C f"^ we obtain \n\f R (t)\ ~ —D\ z t 2 J2 n v n- For distributions with \i < 2 the 
sum v\ is dominated by the largest t> n 's, and thus, In |/r(£)| oc —D\ 2 t 2 f max - Finally, 
we can calculate the distribution function for the strength of the strongest fluctuator 
f max and obtain 

PW = ^ e-^^ . (28) 

fmax 

Most of the weight of this distribution is around w max ~ V- Thus, in a typical sample for 
t < rf x the decay is given by In oc —(D\ z rji) 2 rather than by ([77]) . To understand 

the difference we note that the average decay law ([27]) can also be obtained by averaging 
the realization-dependent — D\ 2 f max £ 2 (valid for t < f max ) over the distribution ([28)1 . 
This average is dominated by rare samples with a fluctuator of strength f max ~ \jt 
rather than by typical samples. 

For longer times, t > 7"^, the integration gives 
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for 1 < fi < 2 

Wn(t)\)F oc lnW 

for fi < 1 



ln(7max/7min) 



(Dx,z¥Y • (29) 



(In \f R {t)\) F cc -c(D AiZ r// 7max )^ 7max t . (30) 
Both results are not self-averaging. 

8.2. Echo signal decay 

For short times, t < J^ x , we find 

(ln\f E (t)\) F oc -c(D XtzV r lma J 1+ » , (31) 

For c^Dx. zV > 7max the echo decay is dominated by this quasi-static contribution; the 
decay takes place on the time scale shorter than the flip time of the fastest fluctuators, 
^/imax- In this regime (c^^Dx^T) > 7 ma x) the result is self averaging since it is 
dominated by fluctuators with D x , z v ~ (cD^ z 7f^^ maiX ) 1/( - 1+ ^ < c 1/>M Dx,zV < D x , z v^ x . 

For longer times, t > the dephasing is due to multiple flips of the fluctuators. 
These times are relevant if c^^Dx^V < 7max- The decay law is given by 

• for 1 < /j, < 2 

(ln|/ E (t)|) F cx-c( J D A ^)'', (32) 

• for fi — 1 

(In \f E {t)\)F ex -c {Dx, z r]t) ln( 7max t) , (33) 

• for fi < 1 

(In \f E {t)\) F OC -ciDx^/lmaxYlmaxt ■ (34) 

All these results are not self-averaging. 



9. Quadratic coupling 

At the optimal working point, the first-order longitudinal coupling D x , z vanishes. Thus, 
to first order, the decay of the coherent oscillations is determined by the relaxation 
processes and for regular power spectra at low frequencies one expects from Eq. (|13|) 
that T 2 = Tx/2. On the other hand, for power spectra which are singular at low 
frequencies the second-order contribution of the longitudinal noise can be comparable 
or even dominate over Ti/2. To evaluate this contribution, one has to calculate 

f2(t) = (e W j 9(f) X\t')dt^j ^ , (35) 

where for the analysis of the free induction decay (Ramsey signal) one sets g(t') = 1, 
while for decay of the echo-signal g(t' < t/2) = —1 and g(t' > t/2) = 1. 
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1/f noise: The free induction decay for the 1/f noise with a high cutoff uj c (the 
highest energy scale in the problem) has been analyzed in Ref. |21j . Depending on the 
time t, the decay is dominated by low- or high-frequency noise, and the decay law can be 
approximated by a product of low-frequency (w < 1/t, quasi-static) and high-frequency 
(u > 1/t) contributions, f2,n(t) = f^nit) • f^ni^)- The contribution of low frequencies 
is given by £D 122 123 

/£*(*) = = • (36) 
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For 1/f noise the variance of the low-frequency fluctuations is o\ = 2Ahi(l/ui X t). This 
contribution dominates at short times t < [(d 2 ujQi/dX 2 ) A/2] \ At longer times, the 
high-frequency contribution 



oo 



« -t J ^ In (l - 27TZ ^ , (37) 

~i/t 

takes over. When t > \{d 2 uo m /d\ 2 ) A/2] 1 (provided uo c > vr (<9 2 w i/<9A 2 ) A) we obtain 
asymptotically ln/^(t) w — /2)(d 2 Uoi/ d\ 2 ) At. Otherwise the quasistatic result (|3*Bjl 
is valid at all relevant times. One can also evaluate the pre-exponential factor in the 
long-time decay. This pre-exponent decays very slowly (algebraically) but differs from 
1 and thus further reduces f2,n(t) |2lj . 

Quasi-static case: In this case, i.e., when the cutoff u c is lower than 1/t for all 
relevant times, the Ramsey decay is simply given by the static contribution (|36p. At 
all relevant times the decay is algebraic and the crossover to the exponential law is not 
observed. More generally, in the static approximation with a distribution P(SX), the 
dephasing law is given by the Fresnel-type integral transform, 

f!%(t) =Jd(X) P(X) exp (i\ ^ X 2 tj , (38) 

which reduces to Eq. (j3T?j) for a Gaussian P(X) oc exp (— X 2 /2a 2 x ). In general, any 
distribution P(X), finite at X = 0, yields a long time decay of /f^ proportional to 

t- 1 ' 2 . 

For fi < 2 the analysis is technically more complicated. In that case the distribution 
of initial conditions P(X ), and equivalently the sum X = J2n v n cr n,z(t = 0), are no 
longer Gaussian-distributed and, in particular, they cannot be characterized by a typical 
width a, due to the divergence of the second moment (v 2 ) of (jlj). The generalized central 
limit theorem tells us that x = X /r) is then distributed according to a Levy distribution 
L^fi^x) and consequently, according to the free induction decay in the quasi-static 
regime is given by 

fUt) = Jdx L, fi (x) exp ^ v 2 t ■ x 2 ^j . (39) 
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For some values of fi explicit expressions are known. An example is the Cauchy 
distribution, L 10 (x) = 1/[k(1 + x 2 )]. Using (jSUjl the free induction decay in the static 
regime, t < 7 max , is then given by 

/ 2 s ^) = e- 4Qi [l-*(v^/i)] • (40) 
Here we introduced the rate 
1#W 2 

a= 2^A^ 2 (41) 

and &(z) = 2n~ 1 / 2 J* dx e~ x denotes the error function. One can expand $(z) in (|%U|) 
to find the asymptotic behavior of for /j = 1: 

, f 1- at) 1/2 for 

J , ■ ( 4 2) 

[ (vrat) 7 for t > c^ 1 

The initial decay for t <C a -1 is thus very fast, but at times t ~ a" 1 the decay crosses 
over to a much slower power law oc 1/y/i. The dephasing time scales as a -1 , but with 
a relatively large prefactor due to the slow algebraic decay. For other values of /i < 2 
the asymptotic behavior of /J 1 ^ has been obtained in Ref. [21]: 

f 1 - C(u) (atf /2 for t^a" 1 

{ D(fi) (at)~ 1/2 for f>a- 1 
where C(/x) and -D(/i) are factors of order 1. 

Let us now discuss the shape of the decay of /Jj qualitatively and comment on their 
validity. The initial decay of /Jj for /i < 2 is singular and thus very fast compared to the 
Gaussian case /x > 2. This initial decay is dominated by strongly coupled fluctuators, 
i.e., by the tail of the distribution (JU). It is, thus, not self-averaging. 

On the other hand, for longer times, t ^> a -1 , the decay goes over to a much 
slower power law. The exponent —1/2 is independent of /i and coincides even with 
the prediction of the Gaussian model (/i > 2). Hence, the 1/y/t decay law appears 
to be universal in the presence of quasi-static noise, independent of the considered 
statistics. For low enough 7 max such that a ^> 7 max the free induction signal decays 
already in the quasi-static regime, t < 7max> an d is thus given by (]4*H|) . Otherwise, 
further analysis characterizing the contribution of fast fluctuators, 7 > t -1 , is needed to 
describe decoherence. 



10. Conclusions 



We have shown that non-Gaussian 1/f noise of ensembles of two-level fluctuators 
frequently leads to non-self-averaging dephasing laws. Non-Gaussian noise arises, for 
instance, when the distribution of coupling strengths between the two-level fluctuators 
and a qubit has a long algebraic tail. In this case, since experiments are performed on 
specific samples, one should study the typical rather than ensemble averaged behavior. 
Interestingly, in certain regimes, e.g., for short-time echo decay, the decay law is self- 
averaging. 
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